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We describe low-lying collective excitations of atomic nuclei with the multi-reference covariant 
density functional theory, and combine them with coupled-channels calculations for heavy-ion fu¬ 
sion reactions at energies around the Coulomb barrier. To this end, we use the calculated transition 
strengths among several collective states as inputs to the coupled-channels calculations. This ap¬ 
proach provides a natural way to describe anharmonic multi-phonon excitations as well as a deviation 
of rotational excitations from a simple rigid rotor. We apply this method to subbarrier fusion re¬ 
actions of ®®Ni-|-®®Ni, ®®Ni-|-®°Ni and '‘°Ca-|-®®Ni systems. We find that the effect of anharmonicity 
tends to smear the fusion barrier distributions, better reproducing the experimental data compared 
to the calculations in the harmonic oscillator limit. 
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I. INTRODUCTION 

In heavy-ion fusion reactions at energies around the 
Coulomb barrier, low-lying collective excitations of the 
colliding nuclei during fusion process plays an important 
role in enhancing fusion cross sections as compared to a 
prediction of a simple potential model dHi. These ef¬ 
fects have usually been taken into account in the coupled- 
channels calculations with a coupling scheme based on a 
harmonic vibrator for spherical nuclei or on a rigid ro¬ 
tor for deformed nuclei ii- With this approach, the 
energy of the first excited state as well as the coupling 
strength from the ground state to the first excited state, 
which can be taken from the experimental data, specify 
all the other excitation energies and coupling strengths 
for higher members in the coupling scheme. Typical ex¬ 
amples which show the subbarrier enhancement of fu¬ 
sion cross sections include the fusion of ^®Ni-|-®°Ni and 
®^Ni-|-®^Ni, for which multi-phonon excitations have been 
shown to play an important role [1, Q- See also Refs. 

[1, 0 for discussions on multi-phonon excitations. Mul¬ 
tiple excitations within the ground state rotational band 
also play an important role in most of fusion reactions 
involved with heavy deformed nuclei ra¬ 
in reality, however, most of atomic nuclei have neither 
a pure harmonic oscillator spectrum nor a pure rigid body 
rotational band. For example, the ®®Ni nucleus, which 
has usually been considered to be a typical vibrational 
nucleus, does not exhibit a level spectrum characteristic 
to the harmonic vibration, that is, the degeneracy of the 
two-phonon triplet is considerably broken. Moreover, a 
recent theoretical calculation based on a multi-reference 
density functional theory with the Skyrme interaction 
also indicates that the B{E2) strengths among the collec¬ 
tive levels in ®®Ni deviate largely from what are expected 
from a simple harmonic oscillator [HI- It is therefore 
of considerable interest to investigate the role of anhar¬ 
monicity, that is, the deviation from the harmonic limit, 


in subbarrier fusion of ^®Ni. 

In Refs. [13)rail the effect of anharmonicity on subbar¬ 
rier fusion of has been investigated using 

the vibrational limit of interacting boson model (IBM). 
See also Ref. ra for an application of this method to 
large-angle quasi-elastic scattering of the ^®0-|-^^^Sm sys¬ 
tem. Although the static quadrupole moment of the first 
excited state could be successfully extracted by analyzing 
the l^h precision experimental data with this approach 
[i3 rai) the application of IBM has several limitations 
for a global study. Firstly, the method is not applicable 
to doubly magic nuclei, as the number of bosons in IBM 
is estimated from the number of valence nucleons out¬ 
side shell closures. Secondly, phenomenological parame¬ 
ters have to be introduced to the model Hamiltonian and 
to the transition operators. It is therefore desirable to 
develop an alternative microscopic approach for nuclear 
collective excitations, which does not rely on the har¬ 
monic limit or the rigid rotor, in order to systematically 
investigate the effect of collective excitations in general 
on subbarrier fusion reactions in a wide mass region. 

In this paper, we employ a beyond-mean-field method 
to describe low-lying collective excitations and combine it 
with the coupled-channels approach to heavy-ion fusion 
reactions. The pure mean-field approximation breaks the 
rotational symmetry and does not yield a spectrum of 
nuclei. This can be overcome by going beyond the mean- 
field approximation, in particular, by carrying out the 
angular momentum projection. One can also take into 
account the quantum fluctuation of the mean-field wave 
function by superposing many Slater determinants with 
the generator coordinate method (GCM). When the pair¬ 
ing correlation is important, the particle number projec¬ 
tion can also be implemented. Such scheme has been 
referred to as a multi-reference density-functional theory 
(MR-DFT), and has rapidly been developed in the nu¬ 
clear structure physics for the past decade [lira- 

The paper is organized as follows. In Sec. II, we briefly 
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review the coupled-channels approach for heavy-ion sub¬ 
barrier fusion reactions. In Sec. Ill, we present the re¬ 
sults of MR-DFT calculations for the ®®Ni and ®°Ni nu¬ 
clei. To this end, we use the covariant density functional 
theory (CDFT), based on the relativistic framework. In 
Sec. IV, we combine the coupled-channels calculations 
with the MR-CDFT approach. We apply this method to 
subbarrier fusion reactions of ®®Ni-|-®®Ni, ®®Ni-|-®°Ni, and 
“^^Ca-l-^^Ni systems, and discuss the role of anharmonic- 
ity of quadrupole vibrations of ^®Ni and ®°Ni. We then 
summarize the paper in Sec. V. 


II. COUPLED-CHANNELS APPROACH TO 
HEAVY-ION FUSION REACTIONS 

Our aim in this paper is to solve coupled-channels 
equations using inputs from the MR-CDFT calcula¬ 
tions. In principle, one could formulate the coupled- 
channels method fully microscopically using the MR- 
CDFT method. In such approach, the internuclear po¬ 
tentials, both for the diagonal and the coupling parts, 
would be constructed by folding an effective nucleon- 
nucleon interaction with calculated density distributions 
and transition densities [I3- It has been known, however, 
that this double folding procedure fails to work for heavy- 
ion subbarrier fusion reactions That is, one obtains a 
surface diffuseness parameter of around a ^ 0.63 fm when 
a double folding potential is fitted with a Woods-Saxon 
function, whereas experimental fusion cross sections sys¬ 
tematica^ require a much larger value, e.g., a ~ 1.0 
fm [13, . An important fact to notice is that the 

double folding method works only in the surface region 
of the potential. For elastic and inelastic scattering, the 
surface region of the potential is ma inly probed and a 
double folding potential is reasonable (22 h^ . In marked 
contrast, fusion reactions involve both the surface and 
the inner regions, where two nuclei appreciably overlap 
with each other. As a consequence, several dynamical 
effects are important in the inner region [ril - lsTjl , and the 
double folding potential looses its validity. 

Another problem of the fully microscopic formulation 
is that the MR-CDFT calculations seldom yield a perfect 
agreement with experimental data for excitation energies 
and transition strengths, even though an overall agree¬ 
ment is often reasonable. On the other hand, in order 
to describe quantitatively heavy-ion fusion reactions at 
energies close to the Coulomb barrier, it is important to 
use reasonable values for excitation energies and transi¬ 
tion strengths. 

In order to avoid these drawbacks of the fully mi¬ 
croscopic approach, in this paper we employ a semi- 
microscopic approach. That is, we use a phenomenolog¬ 
ical Woods-Saxon internuclear potential and adopt the 
experimental value for the coupling strength between the 
ground state and the first excited state. The coupling 
strengths for higher members are not known well in many 
nuclei, and it is for these values that we employ the MR- 


CDFT calculations, after scaling the calculated values 
with the experimental strength for the transition between 
the ground state and the first excited state. The excita¬ 
tion energies are known for most of the collective levels, 
and we simply use the experimental values whenever they 
are available. 

In the coupled-channels approach to subbarrier fusion 
reactions, one expands the total wave function of the sys¬ 
tem in terms of the eigen-functions of the collective states 
in the target nucleus, as 

^ ( 1 ) 

j 


where r is the relative coordinate between the colliding 
nuclei, and J and L are the angular momentum for the 
target state and the angular momentum for the relative 
motion, respectively. Here, for simplicity of notation, we 
have assumed that the projectile nucleus is inert, but an 
extension is straightforward to the case where both the 
projectile and the target nuclei are excited. We have 
also introduced the isocentrifugal approximation Q , and 
have assumed that the angular momentum for the rel¬ 
ative motion does not change by the excitation of the 
target nucleus. Notice that only the Jz = 0 component 
is excited in the target nucleus in the isocentrifugal ap¬ 
proximation. 

Substituting Eq. o to the projected Schrddinger 
equation for the energy E, that is, — E\'^lml) = 

0, where H is the total Hamiltonian, one obtains the 
coupled-channels equations for the radial wave functions 
uj{r) as[i]. 


' 2fi dr^ 


L{L + l)h^ 


^ Vjj> {r)uj> (r) = 0, 


+ Vb(r) - E + ej 


uj{r) 


( 2 ) 


where /i is the reduced mass for the relative motion, Vq (r) 
is the bare potential, and ej is the energy of the target 
state J. Vjj'(r) are the coupling matrix elements given 

by 


Vjj'{r) = (v3jo|V"coup(?',aAo)|</5j'o), (3) 

where I4oup is the coupling potential and aAo is the ex¬ 
citation operator with a multipolarity A. We solve the 
coupled-channels equations by imposing the incoming 
wave boundary condition at r = Tabs inside the Coulomb 
barrier Q, that is, 


uj{r) 


exp (-i kj(r')dr^ (r < r.,.), 

( 4 ) 

- ^Sjjy+\kjr) (r ^ oo), 

( 5 ) 
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where and -ff are the outgoing and the incoming 
Coulomb wave functions, respectively. Sjj^ and are 
the nuclear 5-matrix and the transmission coefficient, re¬ 
spectively, with Jo = 0 being the spin of the ground state 
of the target nucleus. kj{r) is the local wave number 
given by 




whereas kj = kj{r = oo) = — ej)lh?. The fu¬ 

sion cross section crfus is then obtained as 


withPz,(j;) = x;,/|T/j„P- 

As we have mentioned, we employ the Woods-Saxon 
potential for the nuclear part of the bare potential, 
bo^^^(r), that is. 


yW(r) 


1 -b exp[(r - Ro)/a]' 


( 8 ) 


The nuclear coupling potential is obtained by deforming 
the radius Rq to 


i?o —> i?o + Rt (9) 

A* 

where Rt is the radius of the target nucleus. Here, the 
deformation parameter a\fj_ is related to the electric mul¬ 
tipole operator as (3^ . 

3p 

Qxt = ^ZtR^tC^x^, ( 10 ) 

47r 

where Zt is the atomic number of the target nucleus. 
The coupling potential in Eq. ([3|) in the isocentrifugal 
approximation then reads 


For the quadrupole harmonic oscillator with A = 2, this 
definition yields @ 

=72 4^=^^(20201 JO), (13) 

for the coupling between the one phonon state (that is, 
the state) to the two phonon state with the angu¬ 
lar momentum J. Notice that Y^X]j=o 2,4(/5j27^^)^ ~ 

, which has often been employed in the coupled- 
channels calculations with multi-phonon couplings [3, H, 
01. For a rigid rotor, one obtains 

(14) 

with 44^^ = (®®® (3.49) in Ref. @). 

With microscopic nuclear structure calculations, the 
coupling strengths can be estimated as (see Eqs. m 
and (dH), 



2A -b 1 47r 


47r SZtcRj 


(y^JolQxol^fJ’o), 


(15) 


where the quantity (</5jolQAol‘b’J'o) can be evaluated mi¬ 
croscopically using the operator Qxfj, = fi^x^i^i)- 
In coupled-channels calculations for heavy-ion reac¬ 
tions, one sometimes uses a different value of the nuclear 
coupling strength from the Coulomb coupling strength, 
see e.g., Ref. 0. While the Coulomb coupling strength, 
Pc, can be estimated from a measured electric transi¬ 
tion strengths, B{EX), the nuclear coupling strengths, 
/3jVj are taken rather arbitrary. One of the big advan¬ 
tages of the semi-microscopic approach is that the nu¬ 
clear coupling strengths can also be estimated by using 
the isoscalar operator for that is, 

Qxt = E = -^AtR^tC^x^, (16) 

i€p,n 


^OUp(^;Q^Ao) — 


1 -b exp[(r - Ro - \/^^RTaxo) 


2X + 1 


ZpZj^c' 


/2A + 1 


T’A+I 


47r 


-avo 


where we have also included the Coulomb coupling po¬ 
tential. The last term is to avoid the double counting in 
Eq. @. 

The matrix elements of Woup, Eq. ([3|), can be eval¬ 
uated with the method employed in the computer code 
CCFULL ii- To this end, one needs the matrix elements 

of the operator In the following, we define 

the coupling strengths for the coupled-channels calcula¬ 
tions as 


whereas the Coulomb coupling strengths are related to 
the EX operator, 

QIE = 6^4^77) = ^ZrRlax^,. (17) 

i€p 


From these equations, one obtains 


Pn _ ( Mn \ 

Pc V ) ' 


(18) 


where M^/Mp is the neutron-to-proton ratio for a tran¬ 
sition given by 


Mn ^ {‘PJM\J2i(^ur^Yxt,{ri)\(pj'M') 
Mp {‘PJM\T,^ep^^Yxp{ri)\(pJ,M')' 


Notice that for a pure isoscalar transition, Mn/AIp is 
reduced to Nt/Zt, and the nuclear and the Coulomb 
coupling strengths are identical to each other, that is, 

Pn = Pc- 
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III. MULTI-REFERENCE COVARIANT 
DENSITY FUNCTIONAL CALCULATION FOR 
®®Ni AND ®°Ni 

Let us now carry out the MR-CDFT calculations for 
the ®®Ni and nuclei and obtain inputs for the 

coupled-channels calculations. These nuclei have been 
studied recently in Ref. using the non-relativistic 

MR-DFT method. Here we repeat similar calculations 
with the relativistic framework, in order to check the de¬ 
pendence of the conclusions on a choice of energy density 
functional. 

In the MR-CDFT [l^, l33-l35l|. the wave function for 
nuclear low-lying states is constructed as a superposition 
of projected mean-field states corresponding to different 
deformations parameter /?, 

\aJM-,NZ) = (20) 

/3 

where a = 1, 2,... distinguishes different collective states 
with the same angular momentum J. Here, |‘i>(/3)) are 
the mean-held states generated by the deformation con¬ 
strained relativistic mean-held (RMF) method. The pair¬ 
ing correlation is taken into account in the BCS approx¬ 
imation. For simplicity, we have assumed the axial and 
rehection symmetries for the mean-held states. Thus, the 
K quantum number (that is, the projection of angular 
momentum onto the z-axis) is zero, and the index K has 
been dropped. We implement both the particle-number 
and the angular-momentum projections. The projection 
operators , and P^ in Eq. (1^01) project the mean- 

held states onto the states with good neutron and proton 
numbers N, Z as well as a good angular momentum J. 

The weight coefficients /^(/3) in Eq. (EUl) are deter¬ 
mined by solving the Hill-Wheeler-Griffin (HWG) equa¬ 
tions [3^ . 

E [^oo(^, ^') - (/5, P')] fi (/3') = 0, (21) 

P' 

where Xo(/3,/3') = ($(/3)|Po{,XP^| $(/?')) and 

= (<I)(/3)|PPg{)P^P^|$(/3')) are the norm and 
the energy kernels, respectively. The prescription of 
mixed density is adopted for the energy kernel. The de¬ 
tailed expressions for the kernels can be found in Refs. 
[H, The solution of the HWG equation pro¬ 

vides the energy levels and the information on the ma¬ 
trix elements ((/?jo|Qao|v^J'o) in Eq.([Tni) with \(Pjm) = 
\aJM;NZ), which are needed for the coupled-channels 
calculations for fusion cross sections. Notice that in most 
of cases it is hard to determine experimentally the sign 
of this matrix element even if information is sometimes 
available on its absolute value from the electric transition 
probabilities Q (see also Ref. 0)- An advantage of the 
present semi-microscopic approach is that the matrix ele¬ 
ments can be estimated theoretically, including their sign 
as well. 
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FIG. 1: (Color online) (a) The total energy for the angular- 
momentum-projected states with J = 0, 2, 4 and 6 for ®®Ni 
as a function of the intrinsic mass quadrupole deformation 
P of the mean-field states. The particle number projection 
has also been implemented. Those curves are obtained with 
the projected CDFT method with PC-PKl interaction [3^ . 
The low-lying collective levels, obtained with the configura¬ 
tion mixing calculation, are also plotted at their average de¬ 
formation p. (b) The collective wave functions given by Eq. 
(1231) for the states indicated in the figure as a function of 
deformation parameter p. 

Figure [T] (a) shows the energy curves of particle- 
number and angular-momentum-projected states with 
J = 0,2,4 and 6 for ^SNi as a function of the intrinsic 
mass quadrupole deformation P of the mean-field states. 
The PG-PKl parameter set is used for the nucleon- 
nucleon interaction. The hgure also shows the low-lying 
collective states obtained by mixing the symmetry con¬ 
served states. These states are plotted at their average 
deformation /3 defined as 

P{Jo.)=Y,P\9im\ ( 22 ) 

p 

where the collective wave functions g^iP) are related to 
the weight functions in Eq. o as 

giiP) = E[-Xo]'/"(/3,/3')/i(/3')- (23) 

P' 

These collective wave functions are plotted in the lower 
panel of Fig. [TJ 

In order to facilitate the discussion on the properties 
of these states, we collect them and make a comparison 
with the experimental spectrum in Fig. [21 The results 
for the ®°Ni are also shown in Fig. [31 One can see 
that the main feature of the energy spectrum and the E2 
transition strength from 2^ to Oj'’ are reproduced rather 
well. These results are qualitatively similar to the re¬ 
sults of the previous MR-DFT calculations with the non- 














5 


8 


6 

> 

0 ) 

S 4 

X 

LU 

2 


0 


(e^fm'') 
6 "_ 

4^_ 

126(8) 

Exp. 




4 ; 


T °3 .3 

270 206 8'2 229 

2 :^-^ . 

150 




PC-PK1 


FIG. 2: (Color online) Comparison of the experimental and 
the calculated low-lying energy spectra of ®®Ni. The experi¬ 
mental data are taken from Refs. , while the calculated 

spectrum is obtained with the PC-PKl force. The E2 tran¬ 
sition strengths are given in units of e^fm'^. 
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FIG. 3: (Color online) Same as Fig. [H but for ®°Ni. 


relativistic Skyrme SLy4 interaction [m, but quantita¬ 
tively the present calculations with the relativistic DFT 
reproduce the experimental data slightly better. For in¬ 
stance, the B{E2) value from the 2^ to the states in 
^®Ni was 248 fm^ in the previous calculation[l2, while 
it is 150 fm"^ in the present calculation, that is much 
closer to the experimental value of 126(8) fm"^ [^.[3^. 

It is interesting to notice that the overall pattern of 
B{E2) values is quite different from what would be ex¬ 
pected for a harmonic vibrator, even though the excita¬ 
tion energies of the 4^, 2 ^, and Oj states are about twice 
the energy of the 2^ state. In particular, the E2 transi¬ 
tion from the 0 ^ to the 2f states is much smaller than 
that from the 4f and the 2t states to the 2f state, which 
is similar to Cd isotopes [i^, Hlj. Instead, the O 2 state 
has a strong transition from the 2 ^ state, which clearly 
indicates that the 0 ^ state is not a member of the two- 
phonon triplet. Notice that the collective wave function 
for each of the 4f and 2^ states have a similar structure 
to one another (see Fig. [IKb)). On the other hand, the 
0 ^ state has a considerably different wave function from 


those states, being dominated by the mean-field configu¬ 
rations around /3 = 0.3. 

Compared to the Oj state, the E2 transition strength 
from the to the 2f states is much larger and is com¬ 
parable to that from the 4^ and the 2^ states to the 2f 
state. This fact makes the O 3 state a better candidate 
for a member of the two-phonon triplets, even though 
the excitation energy is a little bit large and its collective 
wave function is much different from that for the 4^ and 
the 2^ states. A similar conclusion has been reached also 
with the non-relativistic DFT [lTI |. 

Notice that, in the harmonic oscillator limit, the B{E2) 
value from any of the two-phonon triplet states to the 2f 
state is exactly twice the B{E2) value from the 2f state 
to the ground state [42|. The calculated B{E2) values 
shown in Figs. [5] and [31 together with the strong tran¬ 
sition from the 2 ^ to the 0 ^ states, indicate a presence 
of large anharmonicity in the quadrupole vibrations in 
^®Ni and ®°Ni. That is, the calculated B{E2) values are 
significantly quenched from the values in the harmonic 
limit. This fact also implies that the present approach 
for fusion reactions with MR-CDFT provides a natural 
truncation scheme in the coupled-channels calculations, 
whereas the truncation of the phonon spectrum has to be 
introduced in an ad-hoc way if one employs the harmonic 
oscillator couplings. 

Another clear indication of anharmonicity is a finite 
value of quadrupole moment of the first 2+ state. Ex¬ 
perimentally, the spectroscopic quadrupole moment of 
the first 2+ state has been measured to be Q(2,f) = 
—10 ± 6 efm^ for ®®Ni and Q{2,^) = +3 ± 5 efm^ for 
®°Ni [ 4 ^. The present MR-CDFT calculations yield 
g(2+)=-b7.96 and -blO.4 efm^ for 5®Ni and ®°Ni, re¬ 
spectively. Even though the sign of quadrupole moment 
is opposite to the experimental data for ^®Ni, the MR- 
CDFT calculations predict a similar absolute value of 
quadrupole moment to the experimental value both for 
^®Ni and ®°Ni. Notice that the average deformation for 
the 2^ state is very small for both of these nuclei due to 
a large cancellation between the prolate and the oblate 
components, as shown in Fig. [I] A more careful nu¬ 
merical treatment of the calculations would therefore be 
necessary in order to reproduce the correct sign of the 
quadrupole moment, although it is beyond the scope of 
the present paper. 


IV. FUSION OF Ni ISOTOPES 

In the previous section, we have seen that both ®®Ni 
and ®°Ni do not show a typical behavior of harmonic 
oscillator. Let us now investigate how the deviation of the 
spectrum from the harmonic limit affects the subbarrier 
fusion reactions of Ni isotopes. 

We first consider the fusion reaction of two ®®Ni nuclei. 
Since we assume the axial and reflection symmetries in 
the present MR-CDFT calculations, at this moment we 
are unable to describe both the octupole vibration, 3“, 
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TABLE I: The Coulomb coupling strengths for the quadrupole transitions in ®®Ni estimated with the microscopic MR-CDFT 
calculations with PC-PKl force (see Eq. m)- The values in the parenthesis are the corresponding nuclear coupling strengths 
obtained with Eq. (HI). The radius parameter of ro=1.06 fm is used, and the calculated values are scaled with the Coulomb 
coupling strength for the transition from the ground state to the first 2'^ state, which is estimated with the measured B{E2) 
value. 
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TABLE II: The strengths for the Coulomb coupling in ®®Ni 
in the harmonic oscillator limit. Here, we have assumed that 
the third 0^ state belongs to the two-phonon triplet. 
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and the state in the three-phonon multiplets. We 
therefore consider only the quadrupole two-phonon ex¬ 
citations in each ®®Ni nucleus. Table I summarizes the 
coupling strengths, obtained with the MR-CDFT calcu¬ 
lation discussed in the previous section (see Eq. (fTSl) '). 
We use Rt = 1-06 x 58^/^ fm for the radius of ®®Ni. The 
phase of each of the collective wave functions is chosen so 
that the sign of the off-diagonal components is identical 
to that in the harmonic oscillator limit. The Coulomb 
coupling strength for the transition between the ground 
state and the first 2+ state is estimated to be /3 = 0.223 
using the measured B[E2) value with the same radius pa¬ 
rameter. The calculated values of the coupling strengths 
shown in Table I have been scaled to this value, which 
amounts to multiplying a factor of 0.916 to all the cal¬ 
culated coupling strengths. In Table I, the values in the 
parenthesis are the nuclear coupling strengths, calculated 


with the theoretical value for Mn/Mp based on Eq. (1181) . 
For a comparison, we also show in Table II the Coulomb 
coupling strengths in the harmonic oscillator limit, as¬ 
suming the third O^ state to be a member of the two- 
phonon triplet. As we have mentioned in the previous 
section, the coupling strengths in ®®Ni shown in Table I 
reveal some similarity to the harmonic oscillator. That 
is, the coupling strengths from the one phonon state (the 
2f state) to the two-phonon states (that is, the 2^, 4+, 
and 0;^ states) are close to those in the harmonic limit. 
However, there are also pronounced deviations from the 
harmonic limit as well. That is, the strong couplings are 
present between the 2^ and the Oj states, and also be¬ 
tween the 0^ and the 2^ states. The latter is zero in the 
harmonic limit, and so is the former unless the Oj is a 
member of the three-phonon multiplets. 

Figures m (a) and [3] (b) show the fusion cross section 
o'fus(E) and the fusion barrier distribution D{us{E) de¬ 
fined as 0,113 


Dfus{E) 


d'^{Eafus) 
dE"^ ’ 


(24) 


for the ^®Ni-|-^®Ni reaction, respectively. We use the 
Woods-Saxon potential with V))=170.2 MeV, tq = 1.0 
fm, and a=0.9 fm, where the radius Rq is given as 
Ro = + ^T^)- The fusion barrier distributions 

are obtained with the point difference formula with the 
energy step of AEc.m. = 2 MeV in order to be consis¬ 
tent with the experimental barrier distribution extracted 
in Ref. The dashed line shows the result of the 

coupled-channels calculations including up to the double 
phonon states in the harmonic oscillator limit. All the 
mutual excitations between the projectile and the tar- 
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E (MeV) 

c.m. ^ ' 

FIG. 4: (Color online) The fusion cross sections (the upper 
panel) and the fusion barrier distributions (the lower panel) 
for the ®®Ni+®®Ni system. Here, the fusion barrier distribu¬ 
tion is defined as Dfus(i?) = (Eafus)/dE^. The dashed line 
is the result of the coupled-channels calculations including the 
double quadrupole phonon excitations in each ®®Ni nucleus in 
the harmonic oscillator limit, while the solid line is obtained 
by including the 0)*", 2f, 0^, 2^, and 4f states with the cou¬ 
pling strengths shown in Table I. The dotted line denotes the 
result in the absence of the channel couplings. The exper¬ 
imental data are taken from Ref. [45l| for the fusion cross 
sections and from Ref. for the fusion barrier distribution. 


get nuclei are included. On the other hand, the solid 
line in the figure is obtained with the coupling strengths 
shown in Table I. To this end, we include the 0^, 2^, 

2^, and 4f states in the coupled-channels calculations. 
Again, all the mutual excitation channels are taken into 
account. We use the experimental excitation energies 
for these states, that is, e =0, 1.454, 3.531, 2.775, and 
2.459 MeV, respectively. For a comparison, the figure 
also shows the result of no-coupling limit by the dotted 
line. One can see that the calculations in the harmonic 
limit overestimate fusion cross sections at the two lowest 
energies, while the MR-CDFT calculations underpredict 
fusion cross sections around 95 MeV. On the other hand, 
for the energy dependence of fusion cross sections, shown 
in terms of fusion barrier distribution in the lower panel 
of the figure, the MR-CDFT calculation leads to a minor 
improvement by considerably smearing each peak. 

In order to understand the origins for the smearing in 
fusion barrier distribution, we repeat the same calcula¬ 
tions with three different coupling schemes. In the first 



E (MeV) 

FIG. 5: (Golor online) The fusion barrier distributions for the 
®®Ni-|-®®Ni reaction obtained with several coupling schemes. 
In all the panels, the solid line denotes the result of the MR- 
CDFT calculation shown in Fig. [T] In the top panel, the 
dashed line is obtained by setting the energy of the two- 
phonon triplet to be exactly twice the energy of the first 2"*' 
state. In the middle panel, the dashed line is obtained by set¬ 
ting all the diagonal couplings to be zero, while in the bottom 
panel it is obtained by setting the coupling strength to be zero 
for the transition from the second 2^ state to the third 0’*' 
state. 


scheme, the energy of the two-phonon triplets is set to be 
exactly twice the energy of the first 2+ state. The result 
for this is shown in Fig. [S] (a). In the second scheme, 
all the diagonal couplings are set to be zero, while in 
the third scheme the coupling strength is set to be zero 
between the 2^ and the Oj states. The results of these 
schemes are plotted in Figs. [5] (b) and [5] (c), respec¬ 
tively. One can see that all of these coupling schemes 
lead to a more structured barrier distribution than the 
full MR-CDFT calculations, and thus all of these three 
effects, together with the quenching of the coupling be¬ 
tween the one-phonon and the two-phonon states, coher¬ 
ently contribute to the smearing in the barrier distribu¬ 
tion. Among them, the effect shown in Fig. [5] (c) seems 
to yield the largest effect. 

In connection to Fig. [S] (b). Fig. [5] shows the sen- 



















E (MeV) 

c.m. ^ ' 

FIG. 6: (Color online) The fusion cross sections (the up¬ 
per panel) and the fusion barrier distributions (the lower 
panel) for the ®®Ni-|-®®Ni system obtained with the MR-CDFT 
method. The solid line is the same as that in Fig. [4l while the 
dashed line is obtained by inverting the sign of the quadrupole 
moment of the first 2"*' state. The experimental data are taken 
from Refs. [1, 


sitivity of the fusion cross section and the fusion barrier 
distribution to the sign of quadrupole moment of the first 
2+ state. The solid line is the result with the coupling 
strengths shown in Table I, while the dashed line is ob¬ 
tained by changing the sign of the quadrupole moment of 
the first 2~^ state. One can see that the effect of the sign 
of the quadrupole moment is not large, but is certainly 
not negligible. Therefore, the conclusion in Refs. 
remains the same, that is, the sign of the quadrupole mo¬ 
ment of an excited state can be determined with heavy- 
ion subbarrier fusion reactions, when high precision ex¬ 
perimental data are available. For the ^®Ni-|-®®Ni system 
shown in Fig. [51 the experimental data are reproduced 
slightly better with a negative value of quadrupole mo¬ 
ment of the first 2+ state, which is consistent with the 
experimental observation |43i |. 

Finally let us discuss the effect of the second 0+ state, 
which couples strongly to the second 2+ state (see Table 
I). The dashed line in Fig. [7] is obtained by including the 
second O’*' state in the coupled-channels calculations in 
addition to the two-phonon excitations. On the other 
hand side, the solid line shows the result of the two- 
phonon excitations, which is the same as that in Figs. [H 
inland[71 Despite the strong coupling between the second 
0+ and the second 2+ states, one can see that the main 



E (MeV) 

FIG. 7: (Color online) A comparison of the coupled-channels 
calculations for the ®®Ni-|-®®Ni system with (the dashed line) 
and without (the solid line) the couplings to the second O'*" 
state. The upper and the lower panels show the fusion cross 
sections and the fusion barrier distributions, respectively. The 
experimental data are taken from Refs. [3.14^. 

feature of the barrier distribution remains the same even 
when the second O'*" state is included, although the peak 
structure is further smeared by the second 0+ state. This 
is probably because the second O"*" state is not directly 
coupled to the ground state. 

Let us next consider the ®®Ni-|-®°Ni reaction, for which 
high precision data for fusion cross sections have been 
measured by Stefanini et al. Q. Table III summarizes 
the coupling strengths for ®°Ni obtained with the MR- 
CDFT calculations. The main feature of the coupling 
strengths is similar to that for ®®Ni, although the collec¬ 
tivity is somewhat larger in ®°Ni than in ^®Ni. In par¬ 
ticular, strong couplings between the Oj and 2^ states 
and between the 0;^ and 2^ states remain qualitatively 
the same. In addition, the reorientation term (that is, 
the self-coupling term) for the 2^ is much larger in ®°Ni 
as compared to that in ®®Ni. Fig. [5| shows the results of 
the coupled-channels calculations with the Woods-Saxon 
potential with Vb=154.5 MeV, ro=1.0 fm, and a=0.9 
fm. Since we would like to compare between the har¬ 
monic limit and the MR-CDFT calculations, we consider 
only the two-phonon couplings, excluding the couplings 
to the second O"*" state. Unlike the ^®Ni-|-®®Ni system 
shown in Fig. 21 the fusion cross sections are largely un¬ 
derestimated at energies below the Coulomb barrier by 
this calculation. This is probably due to the elastic two- 
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TABLE III: Same as Table [D but for ®°Ni. 


I 

I' 

Ot 

2^ 

0^ 

2+ 

4^ 


Ot 

0 

0.261 

0 

-0.0101 

0 

0 



(0.288) 


(-0.0170) 



2^ 

0.261 

-0.0616 

0.0478 

-0.148 

0.251 

0.126 


(0.288) 

(-0.0562) 

(0.0390) 

(-0.154) 

(0.267) 

(0.117) 

oj 

0 

0.0478 

0 

-0.390 

0 

0 



(0.0390) 


(-0.367) 




-0.0101 

-0.148 

-0.390 

0.164 

-0.0264 

0.107 


(-0.0170) 

(-0.154) 

(-0.367) 

(0.151) 

(-0.0210) 

(0.114) 


0 

0.251 

0 

-0.0264 

-0.0701 

0 



(0.267) 


(-0.0210) 

(-0.0681) 


ot 

0 

0.126 

0 

0.107 

0 

0 



(0.117) 


(0.114) 





FIG. 8: (Color online) Same as Fig. H but for the ®®Ni+“Ni 
system. The experimental data are taken from Ref. @|. 


neutron transfer process, which is not taken into account 
in the present coupled-channels calculations, as has been 
discussed in Ref. [1]. Despite this, however, one may 
expect that the shape of fusion barrier distribution is 
not much affected by transfer channels, unless a multi¬ 
nucleon transfer process takes place (which is unlikely in 
the ®®Ni-|-®°Ni system). The effect of anharmonicity on 
the shape of fusion barrier distribution is qualitatively 


the same as that in the ^®Ni-|-^®Ni system. That is, the 
anharmonicity largely smears the peak structure in the 
barrier distribution. Even though the agreement with the 
experimental barrier distribution gets worse by includ¬ 
ing the anharmonicity effects, it is interesting to notice 
that the MR-CDFT calculations appear to reproduce the 
more recent (preliminary) data of the barrier distribution 
for the same system |46l |. in a more consistent way than 
the result of the harmonic approximation. 

Lastly, we briefly discuss the ^'^Ca, -I- ®®Ni fusion reac¬ 
tion. Fig. inishows the results of the coupled-channels cal¬ 
culations with the Woods-Saxon potential with Vb=135 
MeV, ro=1.0 fm, and a=0.9 fm. The excitations up to 
the two-phonon states are taken into account in the tar¬ 
get nucleus ^®Ni while the one octupole phonon excita¬ 
tion is included for the projectile nucleus, ^°Ca, in the 
harmonic limit. All the mutual excitations are included 
in the coupled-channels calculation. Since the charge 
product is small for this system, the inclusion of the sec¬ 
ond O’*' state in ®®Ni leads to only a marginal change 
both in the fusion cross sections and in the fusion barrier 
distribution. Both the harmonic limit and MR-CDFT 
calculations reproduce well the experimental fusion cross 
sections [l^. However, in the lower panel of the figure 
one can again see that the anharmonicity effect in ^®Ni 
smears the fusion barrier distribution, leading to a better 
agreement with the experimental fusion barrier distribu¬ 
tion as compared to the results in the harmonic oscillator 
limit. 


V. SUMMARY 

We have proposed the semi-microscopic approach to 
heavy-ion sub-barrier fusion reactions. The basic idea 
of this approach is to combine a multi-reference density 
functional theory (MR-DFT) to a coupled-channels cal- 
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FIG. 9: (Color online) Same as Fig. [H but for the “^^Ca+^^Ni 
system. The experimental data are taken from Ref. [i^ . 

culation. The MR-DFT provides transition strengths 
among collective states without resorting to the har¬ 
monic oscillator model or the rigid rotor model. It also 
provides the relative sign for the transition matrix ele¬ 
ments. These quantities are usually not known exper¬ 
imentally, and the MR-DFT provides important inputs 
for the coupled-channels calculations. The excitation en¬ 
ergies, on the other hand, are often known well, and 
one can simply use the experimental values for them, 
although one could use the results of a MR-DFT calcu¬ 
lation if the excitation energies are not known experi¬ 
mentally. In this paper, instead of carrying out a fully 
microscopic calculation for fusion with a double folding 
potential, we employ the semi-microscopic approach with 
a phenomenological Woods-Saxon potential, since it has 
been known that the double folding procedure does not 
work for subbarrier fusion reactions. The advantages of 
this approach include i) deviations from the harmonic 
limit as well as the rigid rotor limit can be taken into ac¬ 
count, ii) it can therefore be applied also to transitional 
nuclei, which show neither the vibrational nor the rotati- 
nal characters, hi) the sign of the matrix elements can 


be determined, iv) the nuclear coupling strengths can be 
estimated using the neutron-to-proton ratio for a tran¬ 
sition, and v) a natural truncation is introduced in the 
coupling schemes. 

We have applied this approach to the ®®Ni-|-®®Ni, 
58Ni-|-60]\ji, 4°Ca-|-®®Ni fusion reactions. We have 

first discussed the spectrum of ®®Ni and ®°Ni using the 
multi-reference covariant density functional theory. We 
have found that there are both similarities and differences 
between the calculated spectra and those in the harmonic 
limit, even though those nuclei have been considered to 
be typical vibrational nuclei in the literatures. We have 
then discussed the effect of the anharmonicities on the 
fusion cross sections and the fusion barrier distributions. 
We have found that the anharmonicities smear the fusion 
barrier distributions, somewhat improving the agreement 
with the experimental data. 

In this paper, for simplicity we have assumed the axial 
and reflection symmetries for the shape of ^®Ni and ®°Ni. 
This has prevented us from including the octupole and 
the three quadruple phonon excitations in the coupled- 
channels calculations. Although the octupole excitations 
would not affect much the fusion cross sections since the 
excitation energy is large BSi, the three quadrupole 
phonon excitations may perturb the results presented in 
this paper. It would be an interesting future problem, 
even though it may be numerically demanding, to relieve 
the restriction for the symmetries and repeat the calcula¬ 
tions in order to investigate the role of the three-phonon 
excitations. 

It would also be an interesting problem to extend the 
treatment presented in this paper to heavy-ion elastic and 
inelastic scattering. In contrast to fusion reactions, the 
double-folding approach is applicable to these reactions. 
One can therefore develop a fully microscopic approach 
to heavy-ion scattering using the multi-reference density 
functional theory. That approach would also be useful in 
applications to nuclear data, such as those related to the 
problem of nuclear transmutation (49j| . 
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